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SUMMARY 


The dynamic properties at various temperatures of braids impregnated with 
polymer can be measured by using the braid as the suspension of a torsion pendu- 
lum. This report describes the electronic and mechanical design of a torsional 
braid pendulum displacement transducer which is an advance in the state of the 
art. The transducer uses a unique optical design consisting of refracting 
quartz windows used in conjunction with a differential photocell to produce a 
null signal. The release mechanism for initiating free torsional oscillation of 
the pendulum has also been improved. Analysis of the precision and accuracy of 
the transducer indicated that the maximum relative error in measuring torsional 
amplitude was approximately 0.7 percent. A serious problem inherent in all 
instruments which use a torsional suspension was analyzed: misalignment of the 

physical and torsional axes of the torsional member which results in modulation 
of the amplitude of the free oscillation. 


INTRODUCTION 

Torsional braid pendulum displacement transducers have been used to gen- 
erate the dynamic mechanical spectra of polymeric materials over a wide tempera- 
ture range (refs. 1 and 2). Previously, a solid polymeric suspension element 
was part of the torsion pendulum, but because of undesirable effects due mostly 
to temperature changes, the element was changed to a multifilament glass braid 
impregnated with the polymer sample to be studied. Tests have shown that this 
braid makes the undesirable temperature effects insignificant. With the impreg- 
nated braid suspension mounted in a furnace tube so that its temperature can be 
controlled, the pendulum is released from an initial torsional offset which 
initiates free torsional oscillation. Ideally, the pendulum motion should have 
a damped sinusoidal waveform from which the damping coefficient, rigidity, and 
shear modulus of the polymer- impregnated braid can be obtained. Since this 
information is very useful for practical application of various polymers and 
since the braid-suspended pendulum system has decided advantages over other 
methods for obtaining this information, considerable effort has been and is 
still being expended to improve transducers and obtain real-time data analysis 
capability (refs. 3, 4, and 5). The existing pendulum displacement transducer 
used to measure the pendulum oscillation is described in reference 1. A fixed 
and a moving polarizer are attached to the pendulum, and a light beam is passed 
through both into a photomultiplier tube. The output of the tube is a function 
of pendulum displacement. A less bulky, simpler to operate, and more sensitive 
transducer was desired with a better method for initiating the free torsional 
oscillation. The ensuing research effort resulted in a torsional braid pendulum 
transducer which meets and exceeds these criteria and which is considered an 
advance in the state of the art. The transducer uses a unique optical system 
for which a patent has been awarded (ref. 6) . It employs a differential photo- 
cell which was made at Langley Research Center, since it could not be obtained 
commercially. These and other innovations are described in this paper. A 



complete mathematical analysis was also made of the transducer precision and 
accuracy. The transducer, the innovations used in its design, and the results 
of the mathematical analysis are applicable to similar areas of endeavor. 


SYMBOLS 

a = C(jj/(ji3' (see eq. (C20) ) 

b = i/2J 

c damping coefficient of torsional member 

Ci,C 2 /--- constants 
D inside diameter of coil 

d window thickness; also diameter of coil wire; also translation of 0 

to 0 ' 

• constants 
G shear modulus 

I mass moment of inertia 

J torsional polar moment of inertia 

polar areal moment of inertia of the torsional member 
K torsional spring constant or torsional rigidity 

L length of torsional member 

lengths of two parts of compound torsional member 
Sf Lagrangian function 

I length of coil 

I damping constant 

m displacement in figure 10; also mass of pendulum 

N number of layers of wire in coil 

n displacement in figure 10; also refractive index; also number of 

half-periods between peak amplitudes 


2 



fractional overshoot 


radius of curvature of window 
total length of wire in coil 

period of torsional oscillation; also kinetic energy 
component of Ty normal to torsional member 
unperturbed torsion in torsional member 
component of Ty parallel to torsional member 
time 

potential energy 

coordinates of center of gravity of pendulum 
phase angle 

angle between incoming ray and normal to curved surface of window 
angle between plane surface of window and exiting refracted ray 
torsional displacement of pendulum 
initial torsional displacement of pendulum 
0 measured from reference datum other than 9=0 

angle between refracted ray and normal to curved surface of window 
standard deviation of mean value 
angular displacement of window 

angular displacements of torsional member from Y-axis 
natural angular velocity 
damped angular velocity 

dot over a symbol denotes differentiation with respect to time. 





DESCRIPTION OF TRANSDUCER 


The torsional braid pendulum displacement transducer with its related 
hardware and electronics is shown in figure 1. Figure 2 is a cutaway drawing 
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Figure 1.- Torsional braid pendulum displacement transducer 
installed with related hardware and electronics. 


of the mechanical-optical part of the transducer which more clearly shows its 
operation. Light transmitted from the source (14) through a servo-driven win- 
dow (16) is collimated (18), passed through a focusing lens (11) and through a 
second window (24) mounted on the torsional braid pendulum, and brought to 
focus on a differential photocell (26). As the pendulum rotates, the motion of 
the light source image on the photocell produces an error signal which is fed 
into an electronic circuit which drives the galvanometer (13) and its window 
(16) in a direction that restores the light image to its original, or null. 
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Figure 2.- Mechanical-optical part of torsional braid 
pendulum displacement transducer. 


position. The current that drives the galvanometer is proportional to the 
pendulum displacement and is taken as the output of the transducer. 

In figure 2 the light source (14) can be placed so that its filament is 
vertical. With power off, the galvanometer window (16) is at right angles to 
the optical axis. The lens (18) collimates the light source. The focusing 
lens (11) was chosen so that the light beam covers less than half the diameter 
of the second window (24). The light beam is brought to focus on the photo- 
cell (26) by moving the photocell support (3) in or out. The pendulum is made 
long enough so that the window suspended from it is approximately centered with 
the optical axis. Initially, the pendulum window is perpendicular to the opti- 
cal axis, and the small permanent magnet (25) is aligned parallel to the optical 
axis and perpendicular to the field produced by the step coils (19) and (27) . 
When energized, the coils give the pendulum an initial angular offset. Deener- 
gizing the coils allows the pendulum to freely oscillate. When this occurs, the 
focused lamp filament image starts to move across the differential photocell 
(26) producing an error signal. This signal is fed into an electronic circuit 
which drives the galvanometer in the same direction as that of the pendulum 
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window to cause the light spot to remain at nearly the null position on the 
photocell. (The galvanometer moves in the same direction as the pendulxam window 
because of image reversal through the lenses.) With a symmetrical optical sys- 
tem, the galvanometer tracks the pendulum in a one-to-one ratio. A voltage pro- 
portional to the current driving the galvanometer is the transducer output. 

Figure 3 is a schematic of the electrical circuit. The filament supply on 
the lamp is voltage regulated below its rated value to maintain the filament at 
a fixed operating condition and to give it long life. The manual step switch 
energizes the coils which produce the initial torsional displacement of the 
pendulum. However, over a wide range of braid torsional stiffness, a change in 
the permanent magnet may be necessary to obtain the desired angular offset. 

The step B input is for remotely initiating the pendulum oscillation. 
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Figure 4 is the wiring diagram for the electronic circuit. The error sig- 
nal from the photocell is fed into operational amplifier A1 . Circuit gain is 



Figure 4.- Wiring diagram of the electronic circuit. 


controlled by the 10-kf2 potentiometer which is adjusted to give the maximum gain 
with acceptable stability. From there, the signal goes to amplifier A2 which is 
coupled to two transistors. Together these act as a power amplifier supplying 
the galvanometer. The remainder of the circuit, to the right of the 3A line and 
below the dotted line, is essentially concerned with the mechanism for initiat- 
ing the angular offset and with the motor control for maintaining the average 
position of the pendulum. 

The current flowing to the galvanometer is directly proportional to the 
pendulum displacement. A voltage which is obtained linearly from the galvanom- 
eter current is fed into amplifier A3, into amplifier A4, and into the motor 
circuit. Amplifier A3 receives a signal which is proportional to the pendulum 
position and inverts its phase 180° without gain. Amplifier A4 needs only the 
ac component of the signal since it is a differentiator whose output is propor- 
tional to the pendulum angular velocity. Its output is phase shifted 90° and 
provides damping when these signals are summed at the input to amplifier A5 . 

This damping only acts when the step coil is energized and is used to damp out 
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any oscillations in the initial offset. The 10-kf2 potentiometer associated 
with it gives damping control. The 50-kf^ potentiometer biases the direction in 
which the pendulum will be offset and acts as an adjustment to obtain a desired 
initial offset. Amplifier A5 , in conjunction with its transistors, acts as a 
power amplifier supplying the coil. 

A motor system, to which the pendulum is mounted, offers a method of main- 
taining the average pendulum position. This is accomplished as follows. When 
the contacts are open (that is, the relay and step coils are not energized), the 
motor maintains the pendulum at its average null position. The response of the 
motor is so slow that it does not respond to the normal oscillations of the pen- 
dulum. When the coil and relay circuit is open, the field effect transistor 
(FET) acts as a low resistance, and amplifier A6 drives the motor. The motor 
operates at a low power requiring 12 volts and 2 mA. Then a switch closure or 
TTL (transistor- transistor logic) logical zero is applied across pin 6A to the 
ground; this closes the relay which energizes the step coils and biases the FET 
to a large resistance state. This inhibits further motor operation. The step 
coil, however, is now energized and its field acting on the permanent magnet on 
the pendulum causes the pendulum to displace through a specified angle. This 
specified displacement is obtained by summing three signals at the negative 
input of amplifier A5. These signals are obtained from the 50-kJ2 potentiometer, 
from the actual displacement of the galvanometer, and from the angular velocity 
of the galvanometer (controlled through a 10-kf2 potentiometer) which provides 
damping in the step mode of operation. 

The two diodes across the gate and the drain of the FET limit the voltage 
to a maximum of about 0.5 volt. This keeps the FET resistance within a reason- 
ably linear range. 


DESIGN CONSIDERATIONS 

In this section, design considerations are described which led to this 
particular design of the torsional braid pendulum displacement transducer. 
Important requirements for the light source are that the part of the filament 
forming the image on the photocell lie in a plane, that its position not shift, 
that it have long life, and that it have a constant intensity to prevent an 
apparent image displacement because of filament inhomogeneities or filament 
stress distortion. An incandescent filament lamp with a voltage regulator sup- 
plying a stable, underrated lamp voltage fulfilled these needs. 

The next consideration is the lens and window optics. If the system is 
nonsymmetric , that is, unequal image and object distances, the two windows have 
different sensitivities in generating the beam displacement. This may be desir- 
able when the range of swing desired from the pendulum has to be different from 
that of the tracking sensor. However, the combined distance of object and image 
is minimum when they are equal. Shorter focal length lenses, for a given diam- 
eter, collect more light and result in a more compact design. However, as the 
focal length becomes shorter, the refractor window causes more broadening of the 
line source, in the focal plane. This broadening does not cause any problems, 
since the photocell responds to the average intensity value. 
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For much more control and ease of setup, a two-lens system, with one 
acting as a collimator, was chosen. Since the optics preclude placing the 
refractor windows in the collimator space between the lenses, one window is 
placed in the object space and the other in the image space. Their placement 
is not critical provided that they always pass the entire beam. A typical sym- 
metrical setup of an optical system showing image broadening and displacement 
of the image by the refractor window is shown in figure 5. 



Figure 5.- A symmetric, collimated optical system showing 
image formation through quartz windows. 


Since a differential photocell was not available commercially, a standard 
photocell was altered as shown in figure 6. It generates a differential voltage 
at its terminals, which when fed to a differential amplifier results in an out- 
put proportional to the light displacement from the null position on the photo- 
cell. Reference 7 discusses the theory of the differential photocell. 

The coils causing the field that produces the initial offset of the pendu- 
lum are connected in series. The main consideration here is to limit the cur- 
rent they draw from the power supply to a reasonable value. Since the handbooks 
do not supply this coil design information, the formula is derived in appendix A, 
with an example of its use. The coils react with a permanent magnet, located on 
the bottom of the pendulum, to produce the initial offset. It is important that 
this offset be repeatable so that the data reduction is made for the same initial 
amplitude. This may be quite important if the damping coefficient is not con- 
stant. Constant offset amplitude is provided over an extreme range of braid 
rigidities by using three different strength permanent magnets and electrically 
changing the coil field strength. 
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Figure 6.- Drawing of a standard photocell showing alterations 
which make it a differential photocell. 


A high torque galvanometer with good linearity and minimal hysteresis is 
required, and one was obtained that met these criteria. 

Since the magnetic fields associated with the galvanometer and pendulum 
magnet may interact and cause instability, the tank (fig. 2) was made of soft 
iron to act as a magnetic shield. It was also designed so that the gas for 
heating the braid was baffled so as not to disturb the pendulum. A magnetic 
stabilizer was also added to the pendulum to minimize perturbations which were 
discovered after the transducer was designed. After the source of the perturba- 
tions was finally isolated, the cause and effects of this problem were analyzed 
and are discussed in appendix B. 

In a damped dynamic system, the response is a function of the forcing fre- 
quency, the natural frequency of the system, and the damping. Generally, it is 
desired that the transducer follow the system oscillation as exactly as pos- 
sible, that is, in phase, with equal or proportional amplitudes. The forcing 
function is the pendulum motion with frequency of about 1 Hz and an exponen- 
tially modulated amplitude. (Since an exponential input results in a constant 
input-output amplitude ratio and a constant phase shift, the pendulum over short 
time intervals and in the damping range of interest is essentially acting as a 
sine wave input to the tracker transducer system.) A response curve for a sinu- 
soidal forcing function is shown in figure 7. Clearly, the flattest response is 
obtained for a damping coefficient of 0.707. This damping also gives the best 
phase shift characteristics, but this consideration is relatively unimportant 
since the natural frequency of the transducer which is essentially that of the 
galvanometer (about 76 Hz) is much greater than the forcing frequency (1 Hz) and 
very little phase shift occurs. Also, at this forcing frequency, there is con- 
siderable latitude in choosing a damping coefficient that gives a relatively 
flat response. 
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Figure 7.- Dynamic response curves for a damped, second order 
system with a sinusoidal forcing function. 


Hence, the system stability is considered more critical. In a closed loop 
system, high gain is desirable. However, as the gain increases, the effective 
damping decreases and a condition of instability may occur. The circuit sta- 
bility can be investigated by obtaining the open loop gain versus frequency or 
phase angle versus frequency. For a control system with a second order charac- 
teristic, response peaking occurs at the damped natural frequency f^j and a 
90° phase shift occurs at the natural frequency f . From this point the gain 
rolloff approaches -12 dB/octave. If, in the open loop, the gain is greater 
than or equal to 1 when a 180° phase shift occurs, the closed loop circuit is 
unstable. A check of the response of the original circuit (fig. 8) showed that 
it was unstable. In order to maintain high gain and stability, the circuit was 
altered to use velocity feedback to stabilize the circuit. This approach was 
unsuccessful; thus the circuit was modified to resemble a first order system in 
its operational range, which has a -6 dB/octave gain rolloff, does not peak on 
the gain, and theoretically has a 90° phase shift at infinity. The frequency 
response of this system, shown in figure 9, is stable. 
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PRECISION 


Factors affecting precision are, first, those elements which cause a non- 
linear output response for a given input and, second, the errors introduced in 
the data reduction. As was mentioned, the errors from the light source were 
minimized and these errors together with those from the lens system are 
inconsequential . 

To determine what error is introduced by window curvature (nonparallelism 
of the window faces) , a calculation was made to see what curvature would pro- 
duce an error of 1 percent of full scale. Full-scale angular displacement of 
the pendulum is ±15°. For an incidence angle of 15°, the angle of refraction 
through a window with refractive index of 1.5 is 9.936°, as shown in figure 5. 
Thus, the full-scale displacement of the light image caused by window refraction 
is (d tan 15° - d tan 9.936°) where d, the window thickness, is 0.635 cm. 
Therefore, the full-scale displacement is 0.0589 cm, and the 1-percent full- 
scale error is 5.891 x 10”“^ cm. 

Figure 10 shows dimensions of the optical system for which the calculation 
of error due to window curvature was made. All the curvature was assumed on 


WINDOW 



Figure 10.- Diagram for calculating error caused 
by window curvature. 
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one side of the window. Index of refraction of the glass is assumed to be 1.5. 
The angle C of the incoming ray can be described by 


sin C 


0.508 

R 


( 1 ) 


tan ^ 


cos 5 


0.508 



0.2581 


0.2581 

R 


( 2 ) 


(3) 



.508 



.2581 


where R is the radius of curvature. According to the simple relation 
governing refraction, the angle ^ of the ray refracted by the curved window 
surface is 


sin g _ 1 

sin C, 1.5 


(4) 


^ = sin 


( 4 = 5 ^) 


(5) 



■ 


1 


• 2 , 
sin 

2.25 


The displacement n (fig. 10) is then 


n = d tan ^ - d tan ^ = d (tan ^ - tan = 0.635 (tan C “ tan 


( 6 ) 
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Likewise, the angle n of the refracted ray exiting from the plane window 
surface is 

sin (C - g) ^ 1 

sin ri 1.5 



From equations (1) and (5) 


sin (C - C) = sin 


sin 


" (^) - ( W) 


( 8 ) 


The displacement m (fig. 10) is 


m = 1.016 tan p 


(9) 


From equations (6) and (9) , the total displacement is 


m + n = 1.016 tan p + 0.635 (tan C " tan 


( 10 ) 


For a 1-percent full-scale error. 


5.891 X 10 = 1.016 


1.5 sin (^ - g) 


yi - 2. 


25 sin^ (^ - O 


+ 0.635 (tan Z ~ tan ^) 


( 11 ) 
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Substitution of equations (2), (5), and (8) in equation (11) gives 




When Y = 0.508 cm, 



( 15 ) 
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Substitute equation (15) in equation (14) to obtain 


^ = ±-|(r^ - + 0.25s) ^^(±2 \Jr^ - 0.258 ) 

dx 2 ' / \ ^ / 

= 1.969 - 0.258 

Thus , 


tan (90° - C) = = 1.969 UR^ - 0.258 


( 16 ) 


(17) 


cos ^ 
sin ^ 


1.969 



0.258 


(18) 


Substitute equation (3) in equation (18) to obtain 


sin C 



0.258 


R(1.969 



0.258 


1 ^ 0.5079 

1.969R R 


(19) 


Finally, substitute equation (19) in equation (12) to obtain 



For very large values of R, greater than 250 cm, equation (20) becomes 


5.891 X 10 


1 . 524 s sin 


, -1 /0.508\ . -1 / 0.3386 

sxn - srn ( — ^ 


^ 0 ^ ^ 21 ) 


Solving equation (21) for R results in R = 619.8 cm. Thus, a radius of 
curvature of 619.8 cm corresponds to a 1-percent full-scale error. For a 
plano-convex lens 
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where n is the refractive index. For 1/^2 ~ ^1 “ 619.8 cm, the 

focal length is F = 1239.6 cm = 0.081 diopter. 

The curvature of all the windows used in the transducer was checked with 
an optical flat and red laser light with wavelength A of approximately 
700 nm (see sketch) . The applicable equation is 


Y = ± 


\/2 (S.H. 


)R - (S.H. ) 


= ± \ 2 


7.00 X 10" 


R - 


7.00 X 10 


-5 


(23) 



where Y is the distance to the first dark interference ring when the sagittal 
height (S.H.) is A/2. The transducer windows were at least twice as flat as 
those that produced a 1-percent full-scale error. Thus, the radius of curva- 
ture of the transducer windows was greater than 1239.5 cm. Calculating error 
due to window curvature for this value of R (eq. (21)) results in a maximum 
full-scale displacement error of 2.946 x lo~^ cm or 0.5 percent. 


When data reduction is performed as recommended, the nonlinearity of the 
photocell and electronic system output would give insignificant error. The 
error in the peak values of torsional displacement 0 used in the computation 
would conservatively be less than ±0.5 percent. Thus, the total fractional 
standard deviation of the mean in measuring 0, the initial torsional amplitude, 
is 


= \J(0.5 percent)^ + (0.5 percent) = 0.707 percent 


(24) 
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The basic equations used in data reduction are derived in appendix C, an 
the method of obtaining damping coefficient from measurements of three succes 
sive peak amplitudes is described in appendix D. The damping coefficient is 
defined by equation (C21) to be 



where a is calculated according to equation (D15) as 



where 0 q, 6^, and are three successive (all positive or all negative) 

peak amplitudes of the system response curve, as seen in figure 11. 



Figure 11.- Torsional amplitude versus time plot of a damped, 
second order, fully oscillating torsional system. 




The damping error based on the measurement of three successive peak 
amplitudes is now computed as follows : 


dc 


-a^ da 


da 


9 3/2 
(1 + a2) 


o 1/2 
(1 + a2) 


da 

9 3/2 

(1 + a2) 


( 27 ) 


. 8a . 3a 3a 

“ 36q 36^ '^^l 302 ^^2 


27T (0n - 0n ) 


2tt (0„ - 0, ) 


d0n - 


2 

On + 


dc _ da (1 + a^) 
c 


1/2 


, 3/2 

(1 + a^) a 


dc 

c 


2tt (0(-. - 0-I ) 


d0. 


o 

CD 

- 02 

1 — 

CD 

“ 02, 

CD 

O 

- 02 

01 

1 

CD 

da 


+ 

a2) 

o 

CD 

- 02' 


dOj^ + 


''^0 ^1 
i®l “ ®2 


d0. 


0 ®1 \ 2 


O r\ + 

01 - 02 


tl/2 


01 - 02 


d0^ + 


Viii' 

^®1 “ ®2y 


Id0, 


In 


Qq - Qi 

Qi - ^ 2/ \ 

2tt 


1 + 


In 


9q - 01 


2tt 


/ 


(28) 


(29) 


(30) 


(31) 


In equation (31) , it must be remembered that <10 q, d0^, and d02 are 

independent errors even though their amplitudes are dependent on 0q. If all 
the errors are of the same magnitude and sign, then the error dc/c is 0. 

This is not a very reliable estimate of the error. The maximum error would 
occur when d0^^ is opposite in sign to d0Q and d02* The standard deviation 
of the mean, which is independent of sign, is always preferable for giving the 
precision and accuracy of the measurement and is used in the following 
calculation . 
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Let 


be the fractional "overshoot" ; then 



c 



( 32 ) 



(33) 


A practical range of damping coefficients c, taken from reference 1, 
is 0.015 to 0.50. For c = 0.50, the overshoot of alternate peaks is approxi- 
mately 16 percent, and for c = 0.015, the overshoot is 95.4 percent. Calcula- 
tion of maximum error for c = 0.50 with Oq ^9q = 0.707 percent is 


a 

c 

c 


2^/2 (1 + 0.16^ + 0.16^*) 


1/2 


(1 - 0.16“") In 


0.16 


'In 


1 + 


0.16 


2tt 


= 0.212 percent 


(34) 


and the maximum error for c = 0.015 is 


"c ;; 

c 6 




2^^^(1 + 0.954^ + 0.954^) 


1/2 


“ 


(1 - 0.954^) fin — ; 

0.954' 


'In 


1 + 


0.954^ 
2it ) 


- 195 percent 


(35) 



Although a 195-percent relative standard deviation may seem large, it does 
approach infinity as c approaches zero. At that point, O^/c becomes 
indeterminate. For the example given, c = 0.015 ± 1.95(0.015) results in 
0.044 as the upper limit and -0.014 as the lower limit. Since negative damping 
does not exist, the lower limit is zero. Thus the bounds of the damping are 
from zero to 0.044 with a most probable value of 0.015. As damping increases, 
o^/c decreases. If three alternate peaks rather than three successive peaks 
were used, the reader can verify (from eq. (D16) ) that O^^/c would be greater 
than the three-successive-peak method by a factor of 


2 (nn + 1) 




2 n1/2 


(36) 


This factor has a maximum value of 4.0 for n-|^ = 1.0. 

The rigidity or torsional spring constant K is given by equation (C14) 
as 


4tt J 9 2 
K = + c^W J 


= JO) 


'2 


1 + 


(1 - C^) 



(37) 


where J is the torsional moment of inertia, T is the period of oscillation, 
and 0 ) and oj' are the natural and damped angular velocities, and where 
(eqs. (CIO) and (C13) ) 


OJ' 


9 1/2 

(1 - C-^) 00 


and 


00 'T = 2TT 


(38) 


The error dK can readily be obtained as was demonstrated previously for dc, 
by taking 


or 


dK 


9k ^ 3k 
3j 3oo' 


doo' 


,3k 

+ dc 
3c 



(39) 


( 40 ) 
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From equation (37) 


In K = In J + 2 In 0 )' - In (1 - c^) 


(41) 


For the relative error. 


dK _ dJ ^ d(jj' ^ 2c^ dc 

K ~ J (jj' , 2 c 

1 - c 

or 



(42) 


(43) 


Torsional amplitude can be given by 


0 


GJ, 


(44) 


where Tq is the torque in the torsional member of length L, G is shear 
modulus, and is the polar areal moment of inertia of the torsional member. 

Then, since T^/G = K, 



(45) 


For the relative error in shear modulus. 

In G = In J + 2 In co' - In (1 - c^) + In L - In 
d^ ^ ^ ^ 2 dco' ^ 2c^ ^ 

G J CO* _ ^2 c L 


(46) 


(47) 
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DISCUSSION 


The optical system utilizes a newly invented technique for which a patent 
has been obtained (ref. 6) . This consists of the refracting quartz windows 
used in conjunction with the differential photocell for producing a null signal. 
The error contributions from these components can be controlled to insignifi- 
cance, with no undue demands made for their fabrication. This factor together 
with the null signal method of obtaining output, which is probably the most 
advantageous in transducer design, contributes to producing an optimal trans- 
ducer package. An analysis of the electronic circuit did not predict quirks 
that developed during laboratory testing, and numerous additional tests were 
made to optimize the design. 

One of the quirks that appeared, which could not be eliminated by design, 
was a modulation which occurred on the decay curve envelope. A basic analysis 
isolated the cause to any system with torsion in it that had nonstraight-line 
coincidence of the torsional and physical axes of its torsion member. This 
produces a time variation of the torque in the torsion member and a time vary- 
ing coefficient of the angular displacement in the differential equation of 
motion for the system. This time variation causes the modulation of the 
envelope. In a flexible suspension, one of the most likely causes for mis- 
alignment of the torsional and physical axes is rigid end point fastenings. 

These are the most commonly used in practice, and they cause axial misalignment 
of the suspension when it is nontorsionally perturbed. Relatively unrelated 
transducer systems such as those using vibrating wires with clamped ends for 
measuring angular rate have produced odd outputs which can be explained by this 
analysis, since some torsion inadvertently is built into the fixed end point 
fastenings. To minimize or prevent the effects of axial misalignment, care 
should be taken to ensure that the torsional member axis is straight. If end 
pinning devices are attached to it, these should be made colinear with the 
axis. The end fastenings may also be made rigid to torsion but soft to plane 
motion. Finally, unless a perfect end fastening is used (which is unlikely) , 
some form of stabilizer should be used to prevent perturbations from bending 
the torsional axis. 

The detailed error analysis which was performed was for the purpose of 
producing a precise and accurate transducer. However, the results of this 
effort can also be used for selecting the method or methods which give least 
error in the data reduction and for explaining why more scatter occurs over one 
range of the data than over another. Also, until now, it has been believed 
that obtaining G, the shear modulus, of the composite structure of the polymer- 
impregnated braid used for the torsional suspension was not feasible. In equa- 
tion (48) , which gives the relative error of G, all factors contributing to 
the relative error can readily be determined with the exception of dJ^/J^, 
where is the polar areal moment of inertia. Cross sections of the 

polymer- impregnated braid can be obtained by methods similar to those used in 
getting tissue or metallurgical specimens. From these, it is assumed, on a 
statistical basis, that a most probable areal configuration may be obtained. 

It would then be a simple matter to obtain the areal moment of inertia of the 
composite constituents, the relative error of G, and G itself of the com- 
posite braid. By superposition, the G value of the polymer only could then 
be determined. 
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Another factor of interest is the initiation of the free torsional oscil- 
lations of the pendulum. In order to guarantee repeatability of results and 
avoid introducing a forcing function into the initial conditions of the pendu- 
lum release, the initial starting amplitude, as measured from the average pen- 
dulum position, should be kept constant and held sufficiently long that its 
constancy is beyond doubt. Release should then be sharp and should introduce 
no perturbations. The release mechanism in this transducer meets all these 
criteria. 


CONCLUDING REMARKS 

An improved torsional braid pendulum transducer has been designed which is 
considered an advance in the state of the art. The transducer has been 
described and its precision and accuracy have been analyzed. This research 
has resulted in the following findings: 

1. Any torsional system, or system with torsion in it, having misalignment 
of the torsional and physical axes of its torsion member experiences a time 
variation of its torque output when the system is in the dynamic state. This 
results in a time varying torque in the torsion member and a time variable 
coefficient for torsional amplitude in the differential equation of motion for 
the system. The result is a modulation of the normal output of a freely oscil- 
lating, torsional system. To minimize these effects in flexible torsion mem- 
bers, which are normally undesirable, care should be taken to ensure that the 
torsional member axis and its end pinning devices are colinear. A suspension 
system that is rigid to torsional motion but soft to plane motion may be used. 
Some form of stabilizer to prevent perturbations from bending the torsion mem- 
ber axis should be used. 

2. The step input for exciting the torsion pendulum gave a smooth, crisp 
release from the same starting amplitude. 

3. The null method employed in this transducer for obtaining the output 
makes it virtually insensitive to any but torsional motions. 

4. The maximum relative error (due to window curvature) in the torsional 
amplitudes for this system is 0.7 percent. This results in relative errors in 
the measurement of the damping ratio c of 195 percent for c = 0.015 and 
0.212 percent for c = 0.50; this error decreases as the damping ratio 
increases. Since rigidity and shear modulus are dependent on c, more scatter 
in the data should be expected at the lower values of c for all three 
parameters . 

5. To minimize error in data reduction, it is recommended that three suc- 
cessive (all positive or all negative) peaks be read out. 
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6. Making t±ie windows optically flat to one-quarter of the wavelength of 
sodium light would make them a negligible source of error. The main source of 
error would then be the readout of the successive peak amplitudes, and any 
means to improve their accuracy would directly enhance the accuracy and pre- 
cision of the quantities sought. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 25, 1981 
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APPENDIX A 


COIL DESIGN 

The design of the torsional braid pendulum displacement transducer included 
coils which cause the field that produces the initial torsional offset. The 
main consideration in designing the coils was to limit the current they draw 
from the power supply to a reasonable value. The formula for designing the 
coils is derived in this appendix. 

A coil is desired with a given resistance and the dimensions shown in the 
following sketch: 



For the first layer of wire, 


TT (D + d) — = Length of wire in layer 1 


For the second layer. 


(Al) 


TT (D + 3d) — = Length of wire in layer 2 
d 


For the nth layer. 


(A2) 


Tt|^D + (2n - l)c^— = Length of wire in layer n (n = 1,2,3,...) (A3) 

From equation (A3) , an arithmetic progression is obtained whose sum is 

[d + (2N - 


N I t . 

S = -- < TT (D + d) ^ + TT I 


(A4) 


^^{2D + 2Nd) = ^^(D + Nd) 
2d d 


(A5) 
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where N is the number of layers in the coil. This sum is the total length of 
wire of diameter d in the coil. For the coil dimensions in the sketch, 

S = ^ ^ (1.016 + Nd) = 6.811 ^ + 6.704N-^ (A6) 

d d 

The diameter of the coil is also limited to 2.54 cm, that is. 


2.54 

> 

2Nd 

+ D 


2.54 

> 

2Nd 

+ 1.1016 

(A7) 


A coil with resistance of about 55 Q is desired. Try wire with 
d = 0.0254 cm, which has a resistance of 110 per 32 004 cm or 55 per 
16 002 cm. From equation (A6) 

S = N + 6.704n2 = 16 002 (A8) 

0.0254 

6.704n2 + 268. 2N - 16 002 = 0 (A9) 

Equation (A9) can be solved for N, which to a whole number is 33 layers. 

However , 

2Nd + 1.1016 = 2 (33) (0.0254) + 1.106 = 2.69 

which violates equation (A7) . Therefore 0 . 0254-cm-diameter wire is too large. 

Try 32-gage wire with a 0.02032-cm diameter and resistance of 167.3 0, 
per 30 480 cm. To provide coil resistance of 55 Q, 22 layers were required. 

This resulted in a coil diameter of 1.91 cm, which was acceptable. The total 
number of turns is calculated as follows : 

I 2.134 ^ 

d 0.02032 ^ 

105 X 22 = 2310 turns 
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PERTURBATION CALCULATIONS FOR THE TORSION PENDULUM 


This discussion deals with the factors that affect the normal, regular 
decay curve of the damped oscillation of a compound torsion pendulum and the 
output of any system which is put in torsion and whose torsional axis becomes 
bent or deflected. A simple pendulum is first examined, and the final develop- 
ment is for the actual compound pendulum that is used in the transducer design. 
Consider a simple pendulum in torsion as is shown in sketch (a) . If the axis of 


V7777^^777777- 





Sketch (a) 


this pendulum is always a radial extension from the pin point, any torsional 
motion given to the pendulum would not be changed by any perturbation that dis- 
placed it from its statically stable position. This sort of suspension can be 
approached by using a very flexible diaphragm as the pin point for the pendulum 
suspension. This diaphragm is rigid to axial torsion and soft to moments at 
right angles to this torsion. Further, no linear forces acting on the center 
of mass can produce or change the torsion. The more usual situation is that the 
pin point is rigid and a perturbation producing a deflection of the axis changes 
the distribution of the torque from what it was previously. This results in 
modulations and other effects on the oscillation of the torsional pendulum 
(fig. 12). To illustrate the point, one mode of torque distribution is developed 



Figure 12.- Decay curve envelope modulation due to 
perturbations to a torsional pendulum. 
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sufficiently to show that the output is modulated. A second, more general 
torque distribution mode analysis is outlined, but the detailed solution is not 
made in this paper. 


Depending on the torsional spring constant and the pendulum inertia, the 
torque can distribute itself as shown in sketch (b) . The parallel component Tjj 

is Ty cos ip and the normal component Tjj is 
Ty sin ip. For small angles, the parallel component 
is much larger than the normal component, and 
although this development is for the parallel com- 
ponent only, the normal component of torque also 
contributes to the perturbation of the output, since 
it causes rotation of the pendulum. The torsional 
spring constant is K = GJ^/L, and may be assumed 
constant because the effective length L under 
deflection remains essentially constant as does G, 
the shear modulus, and J^, the areal moment of 
inertia. The bending stresses and strains, however, 
do have an effect on these parameters and cause coupling or interactions which 
are not well understood and were ignored in these calculations. 



Assume that the pendulum is in torsion and receives a perturbation which 
moves it off axis to an angle \p . Since the pendulum oscillates because of 
this, \p is a function of time. The parallel component of torsion T|| is 
dependent on ip, and since this is the torsional spring torque in the differ- 
ential equation of motion of the system, it changes the torsional amplitude 0 
from that obtained with the pendulum in axial alignment. This change, of 
course, is an error which should be eliminated. The normal component Tj^ moves 
the pendulum mass normally to the plane of the paper or in a conical motion. 

The change of \p with time prevents a conical path but not the normal rotation. 
This motion also affects the output to produce an error which can be eliminated 
in the same way as the error due to the parallel component is eliminated. 


For free oscillation, the differential equation now controlling 9 in the 
offset axis, considering the parallel torque component only, is 


J0 + £0 + T|| = 0 


(Bl) 


The design of the transducer called for a suspension which results in a 
compound pendulum. The following development for obtaining 0 with perturba- 
tion for this situation is general until the algebra becomes too cumbersome and 
then actual values, taken from the designed pendulum, are used. Assume that a 
torsion Ty has been induced at point B in sketch (c) and a perturbation pro- 
duces and (J). Then, T|| which produces the torsional oscillation is 

Ty cos ip cos (({) - ip) . 
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The conservative Lagrangian is now set up. The only effect of the damping 
force is to reduce 4i and (}) more quickly; it does not alter the argument that 
modulation occurs. Thus, it is omitted in the following determination of ip 
and (j). The position of the center of gravity (point G) in terms of the coor- 
dinate system shown in sketch (c) is 


X = 

sin ip + 

^2 

sin (j) 

(B2) 

Y = L]_ 

cos ip + 

^2 

cos (j) 

(B3) 


The velocity components of the center of gravity at point G are 


• • • 

X = Lj (cos + L 2 (cos (}))cj) (B4) 

y = -j^L]^(sin + L 2 (sin c|))(j)J (B5) 

so that the velocity Vq of point G is 

Vg^ = j^L^(cos + L2(cos + j^L-|^(sin ip)ip 

The kinetic energy of the pendulum is 



- L-5(sin 




(B6) 
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where m is the pendulum mass and Iq is the mass moment of inertia about an 
axis through the center of gravity and perpendicular to the rod. Substituting 
equation (B6) results in 


T = ^ ^1^^ (cos ip)\p + L 2 (cos cj))^ + (sin ip)lp + L 2 (sin 

= + 2L3^L2[cos (cl) - + 






(B7) 


To keep Iq general in terms of m and L 2 , let 
nmL2^ = Iq 

2 

where L, = 25.96 cm. Ip = 4935 g-cm , and m = 38.7 g. Then, 


4935 


38.7 (25.96) ' 


= 0.189 


Substituting Iq = 0.189mL2 in equation (B7) gives 
T = + l.igL^^ij)^ + 2L^L2 jcos (<p - 


(B8) 


The potential energy of the system relative to its equilibrium position is 


V = mg^(Lj^ + L 2 ) “ = mg jl 


- 


cos \ll) + 


L^d 


COS 


4))] 


(B9) 


where g is the acceleration due to gravity. By considering oscillations 
caused by small displacements, cubed and higher order terms are eliminated, and 
since 


, , , (cb^ - 2d)ib + ) 

cos (()) - lb) = 1 - Y 


equations (B8) and (B9) become 


ml'?'') 1' 1 

T = 2di W + I.I 9 L 2 (j) + 21^12# 


(BIO) 
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r 


Consider the Lagrangian, S£ = T - V, with generalized coordinates 4' arid (p 
represented by q^.: 


^ _ _3g_ 

dt 34^ 3q^ 

which may be written as 


(B12) 



dt 


3V 

3qr 


0 


(B13) 


since -3v/3q^ 
to obtain 


3t 

3ij; 

3t 

3(ii 


= 0 and 3T/3q^ = 0. Differentiate equations (BIO) and (Bll) 




(B14) 


^|l.l9 (2)L22(}) + 2Lj^L2li>J = m(l.l9L2^(t) + L^L2’^) 


(B15) 


= ^(2LiI()) = mgLiil; 


(B16) 


dcp 


^(2L2(p) 


mgL24) 


(B17) 


Equations (B14) to (B17) substituted in equation (B13) give the equations of 
motion : 


+ L;LL2^)j + mgLj^i|; = 0 

. 19Ij 2^4’ ^1^2’i' ^ rr>gL24> = 0 


(B18) 

(B19) 
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From equation (B18) , 


+ L 2 $ + gip = 0 


(B20) 


From equation (B19) , 


+ 1.19L2<)) + g(j) = 0 


(B21) 


Equations (B20) and (B21) give 


LiP^ - g 



1.19L2P^ 


g 


0 


(B22) 


where the p's are the angular frequencies of the principal modes of oscilla- 
tion. Thus, from equation (B22), 

0.19L2L2p'^ - (L]_ + 1.19L2)gp^ + g^ = 0 

where p^ is 



(Li 


2„2 


2^ + 1.19L2)g + xj (L]^ + I.I 9 L 2 ) g - 4 (0.19)Lj^L2g 

2 (0.19L^L2) 


+ 1.19L2)g ± \J 


L^^g^ + 1.62Lj^L2g^ + (1.19L2)^g^ 


0. 38L^L2 


(B23) 


From equation (B23), it can be seen that both roots are positive. The periods 
are then 


T 


1 



(B24) 


T 


2 


27T 


P2 


(B25) 
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P 


2 2 

where p-j^ and P 2 are the positive roots of pj^ and P 2 . The angles i 
and 4> are still needed. Therefore, subtract equation (B20) from equation 
to get 


0.19L2cf) + g (<J) - ip) = 0 


'1 r- g 5.26gU) 

(() + 5.26 ()> = 0 

L2 t.2 


Multiply equation (B20) by 1.19 and subtract equation (B21) to get 


Tp + 6.26 2^ - 5.26 <p = 0 


In operator form, equations (B26) and (B27) are 


( d 2 + 5.26 — ]d) - 5.26 ip 0 

-5.26 -3- (j) + (d^ + 6.26 = 0 

Li \ 

The determinant of equations (B28) and (B29) is 


A = 


+ 5.26 -3- -5.26 ~ 


-5.26 


+ 6.26 

L-i 


Let 5.26 r — = M; then 


(e.26 ^ + m)d2 + M A 


A = + [6.26 + m'\d^ + M/6.26 -f- 


= D'^ + 


M 




Let N = 6.26 + M; then 


A = + ND^ + ^ 

^1 


(B21) 

(B26) 

(B27) 

(B28) 

(B29) 


(B30) 

(B31) 
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4 ? 

D + ND“^ + 




-M 


d2 + 6.26-3- 


D** + ND^ + ^j(() = 0 


D"^ + Nd2 + = 


+ M 0 


-5.26 


= 0 


D'^ + ND^ + 


= 0 


For both (p and \p, 


A ^ = 0 

^1 


(B32) 


(B33) 


where r is a characteristic root. . Then 



/_ 5.26 _ 6.26 \ ^ I 39.2 ^ 44.8 ^ 27. ~7 
V ^2 y V Li2 L 1 L 2 ^ 


(B34) 


At this point the following actual pendulum values are substituted into 
equation (B34) : 

L-|^ = 21.84 cm 

L 2 = 25.96 cm 
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= 980 


/ 5.26 

6.26 \ 


39.2 1 44.8 1 27.7 

\ 25.96 

21.84 ) 

N 

(21.84)^ (21.84) (25.96) (25.96)^ 


= 980 


-0.489 ± Q.202 \ 

2 ) 


= -338.6 
= -140.6 
= ±18.41 


(B35) 


(B36) 


r2 = ±11.91 


(B37) 


Equations (B36) and (B37) give the solution to the differential equations. 
Since there is no particular solution, let 


ip = Cj^e ^ 18.4i)t ^ ^^g(11.9i)t c^e ^ 11.9i)t (B38) 


= d,e(18-4i)t ^ (-18.4i)t ^ ^ (11.9i)t ^ (-11.9i)t 


(B39) 


Differentiate equations (B38) and (B39) to obtain 


; (18.4i)t,,o /I ■ \ (-18.4i)t, ... 

ip = Cj^e (18, 4i) + 02^ (-18. 4i) 


+ C3e^^^'®^^^(11.9i) ± c^e^ ^ ^ (-11 . 9i) 


(B40) 


ijj = c^e^^®"'^^^^(-l) (18.4)^ + C2e' (-1) (18 . 4) ' 


+ 036 ^^^‘^^^^(-1) (11.9)^ + c^e[ ^^-^^^^(-1) (11.9) ^ 


(-18.4i)t| 
(-11.9i)t 


(B41) 


(}) = dj^e ^^®-'^^)^(18.4i) + d2e ^ ^ (-18 . 4i ) 


± d3e^^^‘®^^^(11.9i) ± d4e ^ ^ ^ (-11 .9i) 


(B42) 
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^ (18.4)^ + ^ ^ (-1 ) (18 . 4) ^ 


+ (11.9)^ + d^e^ ^ (-1) (11 . 9) ^ 


(B43) 


Substituting equations (B38) , (B39) , and (B43) in equation (B26) gives 


-d^e^^®’^^^^(18.4)^ - d^e^ ^ ^ (18. 4) ^ - d^e ^ ^ (11 . 9) ^ 


. (-11.9i)t 2 , , ,, g r ^(18.4i)t (-18.4i)t ^ , (11.9i)t 

- d e (11.9) + 5.26 :p— d e + d e + d e 

4 J_i2 L -1- ^ 5 


^ (-11.9i)t (18.4i)t (-18.4i)t (11.9i)t (-11 

+ d_e - c e - c_e - c^e - c,e 

4 12 3 4 


•9i)tJ ^ 0 

(B44) 


Substituting equations (B38) , (B39) , and (B41) in equation (B27) gives 


(18.4i)t„„ .,2 (-18.4i)t,^„ ^.2 (11.9i)t,,, ^.2 

-c^e (18.4) - c^e (18.4) - c^e (11.9) 


- + 6.26 ^ 


r 


(18.4i)t ^ (-18.4i)t 

e + c^e 


J 


1.26 ^ d,€ 
L]_ L -L 


(18.4i)t 


. (-18.4i)t ^ . (11.9i)t ^ . (-11.9i)tl 
+ d^e + d^e + d^e = 0 


(B45) 


Equating equations (B44) and (B45) gives 


-d^e^^®'^^^^(18.4)^ - d^e^ ^ ^ (18 . 4) ^ - d^e ^ ^ (11 . 9) ^ 


(-11.9i)t 2 g r (18.4i)t , . (-18.4i)t , . (11.9i)t 

- d.e (11.9) + 5.26 - — d e + d e + d e 

4 ^2 L ^ 


, ^ (-11.9i)t (18.4i)t (-18.4i)t (11.9i)t (-11 

+ d^e - c^e - c^e - c^e - c^e 


.9i) t 


(18.4i)t„„ ^,2 (-18.4i)t„„ .,2 (11.9i)t„, „,2 

= -c^e (18.4) - c^e (18.4) - c^e (11.9) 


(-11.9i)t 2 

c^e (11.9) + 6.26 


Li 


(18.4i)t (-18.4i)t 

e + c^e 


(11.9i)t ^ (-11.9i)t“| .0.9 

:.e + c.e - 5.26 -f- 

3 4 J L 


j r (18. 4i 

rh* 


4i)t 


^ . (-18.4i)t ^ . 

+ d^e + d^e 


(11.9i)t ^ . (-11.9i)t 

+ d^e 


(B46) 
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From equation (B46 ) , 

j^-(18.4)^d^ + 5.26 ^ - 5.26 ^ + (18.4)^c^ 

g g j ”1 (i8.4i) t 

: — c. + 5.26 - — <3.. I e 

^ \ J 


- 6 . 26 r — 


= 0 


-338. 6d^ + 198. 6d^ + 236. Od^ - 198. 6c^ + 338. 6c^ - 280. 9c^ = C 


140.9c 


^1 = 


— = 1.47ct 
96 1 


From equation (B46) , 


[- 


(18.4)^d^ + 5.26 d_ - 5.26 + (18.4)^c, 

2 2 ^ 


- 6.26 -2- c„ + 5.26 -2- d 

Li 2 2 


(-18.4i)t _ - 
e =0 


-338. 6d^ + 198. 6d^ + 236. Od^ - 198. 60^ + 338. ec^ - 280. 90^ = ' 




From equation (B46) 


-(11.9)^d^ + 5.26 ^2 ~ h' ^3 (H-9) 


6.26 -2- + 5.26 ~ d, 

j J_i^ .j 


“ (11.9i)t 

3j® 


= 0 


-141. 6d^ + 198. 6d^ + 236. Od^ - 198. 60^ + 141. 6c^ - 280. 




(B47) 


(B48) 


0 

(B49) 
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From equation (B46) , 


-(11.9)^d + 5.26 ^ - 5.26 c + (11.9)^c 

rr "^2 ^ "^2 ^ ^ 


6.26 


— c + 5.26 e =0 

.1 4 4j 


-141. 6d^ + 198. 6d, + 236. Od^ - 198. 6c^ + 141.6c. - 280.9c. = 0 

4 4 4 4 4 4 


^^4 = 


(B50) 


From equations (B47) to (B50) , only four arbitrary constants remain; these 
inserted in equation (B39) give 

, , (18.4i)t ^ (-18.4i)t ^ (11.9i)t 

cp = 1.47c^e + 1.470^6 + l.lSc^e 


+ l.lSc^e 
4 


(-11.9i)t 


(B51) 


The c's are now evaluated. 

At t = 0, an initial offset is assumed as shown in sketch (d) 



l|; = (}) 


(B52) 


Sketch (d) 


(Note that the most general initial condition should be ip^ ^ (j)^ which would 
make the derivation unwieldy. The assumed condition frequently occurs.) Sub- 
stituting these values into equations (B38) and (B51) results in 


= c^ + C 2 + + c_j (B53) 

(j)^ = 1.47c^ + 1 . 47 c 2 + I.I 5 C 3 + l.lSc^ = ip^ (B54) 
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Differentiating equations (B38) and (B51) results in 


,1, /no (18.4i)t (-18.4i)t 

ip = (18.4i)c^e - (18.4i)c2e 

(11.9i)t (-11.9i)t 

+ (11.9i)c^e - (11.9i)c_^e 


1 /l-i/no /I ■ \ (18.4i)t ^ yi • V (-18.4i)t 

1.47(18.4i)c^e - 1 . 47 (18 . 4i ) c^e 

irm Q-n (11.9i)t nc/m o-\ (-ll-9i)t 

+ 1 . 15 (11 . 9x) c^e - 1 . 15 (11 -9i ) c^e 


At t = 0, 4 j = 4) = 0, so that equations (B55) and (B56) become 


18. 4c^ - 18. 40^ + 11. 9c^ - ll-9c^ = 0 


27.05c^ - 27.050^ + 13.690^ - 13.69c^ = 0 


From equations (B53), (B54) , (B57), and (B58) , 


A = 


1 — t 

1 — 1 

1 

1 — 1 

1.47 

1.47 

1.15 

1.15 

18.4 

-18.4 

11.9 

-11.9 

27.05 

-27.05 

13.69 

-13.69 


= -89.3 




o 




I'o 

1.47 

1.15 

1.15 

0 

-18.4 

11.9 

-11.9 

0 

-27.05 

13.69 

-13.69 


^1 = 


-89.3 


-0.235i|J^ 


(B55) 

(B56) 

(B57) 

(B58) 

(B59) 


(B60) 
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C2 


^3 


^=4 


1 


1 

1 

1.47 


1.15 

1.15 

18.4 

0 

11.9 

-11.9 

27.05 

0 

13.69 

-13.69 


-89 

.3 


1 

1 

^o 

1 

1.47 

1.47 

^o 

1.15 

18.4 

-18.4 

0 

-11.9 

27.05 

-27.05 

0 

-13.69 


-89.3 


1 

1 

1 


1.47 

1.47 

1.15 


18.4 

-18.4 

11.9 

0 

27.05 

-27.05 

13.69 

0 


-89.3 


-0.235i/jg 


0.737i|;^ 


0.737i|Jq 


(B61) 


(B62) 


(B63) 


Substituting equations (B60) to (B63) into equations (B38) and (B51) results in 


ip = - 0.235i|;^e^ 18.4i)t ^ (11 . 9i ) t 


+ 0.7371/J e 
o 


(-11.9i)t 


(B64) 


Since 


= 1.47 (-0.235i|;^)e + 1 .47 (-0 . 235i}J^) e 


+ 1.15 (0.737i|jQ)e + 1 . 15 (0 . 737i|;^) e 


(-18.4i)t 

(-11.9i)t 


cos t 


it , -it 
e + e 

2 


(B65) 


(B66) 
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equations (B64) and (B65) become 

^ = -0.470tjj^ cos (18. 4t) + 1.474 i|Jq cos (11. 9t) 
({) = -0.691iJj^ cos (18. 4t) + 1.6951/;^ cos (11. 9t) 
The output torsion is 


T|| = Ty cos ijj cos (cp - ip) 

Equations (B67) and (B68) substituted in equation (B69) give 

T|| = Ty cos ^-0.470t|;^ cos (18. 4t) + 1.4744^^ cos (11. 9t^ 
X cos ^- 0 . 2214 Jq cos (18. 4t) + 0.221ip^ cos (11. 9t^ 

where 


= K0 

Equation (Bl) , the differential equation for the perturbed condition, thus 
becomes 

j’e + £0 + K0 l^cos ^-0.4704;^ cos (18. 4t) + 1.4744;^ cos (11.9t)J 


X cos 


^- 0 . 2214 ^^ cos (18. 4t) + 0 . 2214 ;^ cos (11.9t)^J = 


Let 


cos ^ 0 . 4704 ^^ cos (18. 4t) + 1 . 4744 ^^ cos (11. 9t)^ = cos f (t) 
cos |^0.2214»q cos (18. 4t) + 0.2214|q cos (11.9t)J = cos g(t) 
Thus, equation (B72) , with the notation of the constants changed, becomes 
6 + 2cto0 + 0)^0 cos f (t) cos g(t) = 0 
where c = £/2Joj and co = \J K/J. 


(B67) 

(B68) 

(B69) 

(B70) 

(B71) 

(B72) 

(B73) 

(B74) 

(B75) 
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The modulation frequency is slow relative to the oscillating frequency; 
therefore, under these conditions, the Wentzel, Kramers, Brillouin (WKB) method 
of solution (from ref. 8) can be applied to solve equation (B75) . Assume that 


0 = A(t) 


i6 (t) 
e 


(B76) 


Then 


0 = Aige^'^ + Ae^^ 


(B77) 


'0 = Ai6(iB)e^^ + Ae^^iB + Ai3e^^ + Aipe^^ + Ae^^ 


(B78) 


Substituting equations (B76) to (B78) into equation (B75) gives 


A + i (a'b + 2a 3 - 2icwA) + a[w^ cos f (t) cos g (t) - + 2icco6^ 


0 


(B79) 

The expression in brackets in equation (B79) is arbitrarily set equal to zero: 


B^ - 2ico)B ~ cos f (t) cos g(t) = 0 


(B80) 


Let 


F (t) 


2 

CO cos 


f (t) 


cos g (t) 


a = 2cco 


(B81) 


With equations (B81) substituted in equation (B80) , 


B^ - iaB - F(t) = 0 


(B82) 


B = 


ia 

2 


+ 



+ F(t) 


(B83) 


Since A is slowly varying, A still slower, and A ~ 0, the following is 
obtained from equation (B79) : 


a’b + A(2B - ia) = 0 


(B84) 
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a6 + 2a(b - 


= 0 


3 


e ^ 


+ 0 
la A 


+ 2 In A =? In c. 


Integrating equation (B86) gives 

i„(b - f ) 

- t) -5 

S = i<=s(^ - f) 


■ 1/2 


and substituting for 3 from equation (B83) gives 


r 

"2 


“2 

A = ±a^i± 

- ^ + F(t)_ 

j = ±05(11)^/^ 

- ^ + F{t^ 


-1/4 


Integrating equation (B83) gives 


B = + 

2 


13 




1 1/2 


+ F(t) 


dt 


1 1/2 


+ F(t) 


dt 


Substituting equations (B90) and (B91) in equation (B76) gives 


0 = iCgCil) 


1/2 


+ F(t) 


-1/4 


exp 


{-f *‘/[- 


F(t) 


1/2 


dt 


Equation (B81) substituted in equation (B92) finally gives 

-1/4 


0 = ±05(11) 


-I- 0 )^ cos £(t) cos g 


(tO 


exp 


l^-ccot ± i /[- (cw)^ + 10^ cos f (t) cos g(t)J 


1/2 


dt 


(B85) 

(B86) 

(B87) 

(B88) 

(B89) 

(B90) 

(B91) 

(B92) 

(B93) 
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Note that the integral in equation (B93) may be written as 
±i(jO J' [j-c + cos f (t) cos g(t)J^^^ dt. 

Refer to equations (B73) and (B74) and let 


= -0.470iJj^ cos (18. 4t) 
a 2 = 1.474i|Jq cos (11. 9t) 

= -0.221\|;q cos (18. 4t) 

02 = 0.221iJj^ cos (11. 9t) 
and let Y = 18. 4t and £ = 11. 9t. Then 


cos f (t) cos g(t) = cos (a^^ + a 2 ) cos (3-[^ + B 2 ) (B94) 


and by series expansion. 


cos f (t) cos g (t) = 1 - 2a^a2 + + 3 p^ + 23^62 + 

+ + 2a2^3i32 + 2a3^a2Bi^ 

+ 4a^tt2B]^B2 + 2a^a2B2^ 

+ a2^Bj^^ + 2a2^B3^62 + «2^B2^) + . . . (B95) 

since is a small angle, the small angle approximation can be made and all 

terms higher than second order can be eliminated. This finally gives 
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2 2 2 

cos f (t) cos g (t) = 1 - 0.0675ijj^ cos 2y - 0.0675i/j^ + 0.742iJ;^ cos y cos e 

- 0.555i)j^^ cos 2e - 0.5551);^^ + 0 . 000338ijj^'^ cos 4y 

+ 0.0003384^^^ + 0.00135t|;^'^ cos 2y + 0 . 00135iJ;^'^ 

. 4 4 

- 0-0700Tj;^ cos 2£ cos £ cos Y “ O.OVOOiJJ^ cos £ cos y 

4 4 

+ 0.0158l(J^ cos 2y cos 2£ + O.OlSSl/^^ cos 2y 

+ 0.0158l|j^^ cos 2e + 0.0158l|J^‘^ + 0.00331ljj^^ cos 4e 

+ 0.0033HI; + 0.0132l|; cos 2£ + 0.00663l|^ (B96) 

o o o 


By gathering like cosine terms and taking TpQ 3.s a small angle, it can readily 
be seen that all terms containing can be dropped with little error. 

Then, 


cos f (t) cos g(t) = 1 - 0.622i^^^ - 0.06754^^^ cos 2 y 


- 0.5554Jq^ cos 2 e + 0.7424^^^ cos y cos £ 


(B97) 


and the integral in equation (B93) is 


/[- 


±ico I -c + 1 - 0.6224; - (0.0675 cos 2 y + 0.555 cos 2 e 


211/2 

- 0.742 cos y cos e)4^Q dt 


Let C 0 = -c + 1 - 0.6224;^^ ; then the integral becomes 


/* r 2~1^/^ 

±ico I <^0 “ (0.0675 cos 2 y + 0.555 cos 2 e - 0.742 cos y cos E)4 ^q dt 
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Expanding the integrand by the binomial expansion gives 


[=e- 


(0.0675 cos 2y + 0.555 cos 2Z - 0.742 cos y cos e)ljj 


O^] 


1/2 


1/2 - 1/2 

- - l/2Cg (0.0675 cos 2y + 0.555 cos 2£ 


^ > , 2 , l/2(-l/2) -3/2 

- 0.742 cos Y cos ^ 


Co [(0.0675 cos 2y 


—I 2 A 

+ 0.555 cos 2£ - 0.742 cos y cos £)j ip^ + . . . 


(B98) 


Again taking the small angle approximation for small and discarding all 

terms higher than second order leaves the integral 


±ico 


/ 


. -(0.0675 cos 2y + 0.555 cos 2Z - 0.742 cos y cos e) 

^ 2 c 1/2 


dt 


Substituting y = 18. 4t and £ = 11. 9t gives 


(B99) 


±io) 


/ 


1/2 




2 c 1/2 

ZC-g 


[o.0675 cos 2(18. 4t) + 0.555 cos 2(11. 9t) 


- 0.742 cos (18. 4t) cos (11.9t)[j \ dt 


(BlOO) 


tiwjc 1/^t 

r 


2c, 


1/2 


0.0675 


sin 2(18. 4t) 


2 (18.4) 


+ 0.555 


sin 2(11. 9t) 
2 (11.9) 


0.742 


sin (18.4 - 11.9)t ^ sin (18.4 + 11.9)t 


2(18.4 - 11.9) 2(18.4 + 11.9) 


+ c. 


(BlOl) 


±ito^[l - c - 0 . 622 i|;^ 2 ^ 1 / 2 ^ _ 


0.0675 


2^1 - c - 0 . 622 i|Jq 2 ^ 


1/2 18.4 


sin 2(18 


.4t)~ 


0.555' 

sin 2 (11 . 9t) 

11.9 

2 


- 0.742 


sin (6.50t) ^ si n ( 30. 3t) 


13.0 


60.6 


+ c- 


(B102) 
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Substituting the above for the integral expression in equation (B93) gives 


0 = ±c 


(±1) + (jS^ cos f (t) cos g(t)J 


-1/4 


exp 


-co)t ± ioj 


7 

[l - c - 0.622^^^ 1 t - 




2 1 - c - 0.622\ii^ 


1^0.00 367 2(18.^)^ ^ 0.0466 

0 (6.50t) ^ sin (30. 3t) 


1/2 


sin 2(11. 9t) 


I 


13.0 


60 . 6 


+ c- 


(B103) 


Since is a small angle, cos f (t) cos g(t) ~ 1 (see eq. (B97)) and 

the integrand remains real until the damping almost becomes critical and 0 of 
equation (B93) is no longer valid since A / 0 in equation (B79) . Normally, 
most torsional operation of this kind is sufficiently far away from critical 
damping operation to make equation (B93) quite valid. However, if the applica- 
tion is such that operation is close to critical, then there are computer pro- 
grams which will solve equation (B75) by substituting 


Z 


= 0 


(B104) 


into equation (B75) and getting 


Z + 2c(jjZ + w‘^0 cos f (t) cos g(t) = 0 


(BIOS) 


Equations (B104) and (BIOS) are sufficient for 
the computer format. 

The following is an outline of the procedure 
for obtaining a more general solution to the above 
problem. Since the above solution demonstrates 
that the perturbation does indeed produce a modu- 
lation on the output, it is sufficient for this 
paper. The more general solution is left for 
future development. 

The angles ijj , (p , Yr snd shown in 

sketch (e) are the same or similar in sense to 
those used in spherical coordinates, and 0 is 
the axial torsional angle. By convention, the 
Z-axis as shown should be negative but it is 
assumed positive. 
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The coordinates of the center of gravity are 
X = sin ip cos Y + L 2 sin (() cos 

y = sin ip sin y + L 2 sin (p sin 

z = cos ip + L 2 cos (j) 

The velocities are 

• • • 

X = -L^ (sin ip) (sin y)Y + (cos y) (cos ip)ip 

- (sin (j)) (sin (cos Q) (cos (p)(p 

y = (sin ip) (cos y)Y "*■ (sin y) (cos \p)^ 

+ (sin (j)) (cos (sin f2) (cos (|))4) 

z = -L^(sin ip)ip - L 2 (sin 

The kinetic and potential energy become 

mx^ , my^ , mz^ , ^(j)^ , ^0® 

2 2 2 2 2 2 

V = mgj^L^d - cos i|^) + L^d - cos cf))^ + ^(Q^ + y2 + 
The Lagrangian then is 


(B106) 

(B107) 

(BIOS) 


(B109) 

(BllO) 

(Bill) 

(B112) 

(B113) 


a = 


m p • • 

T - V = — -L^(sin \p) (sin y)y + (cos y) (cos rp)ip 

• d ^ m I” 

- L 2 (sin ({)) (sin + L 2 (cos fi) (cos + — IL^ (sin i|;) (cos y)y 


+ (sin y) (cos ip)ip + L 2 (sin cf)) (cos (sin J2) (sin <J))(^ 


•~l 2 


+ ^jjL^(sin lp)ip - (sin c())<0 


•n2 l 06 ^ 

- — - — + — — + — - — 


- mgL^ (1 - cos i/j) - mgL^ d - cos (J)) - ^(9^ + y^ + ^2^) (B114) 
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where 



and where interactions and other effects might slightly alter 0, the output 
amplitude due to axial torsion/ but most of the axial torque is assumed trans- 
mitted at all angles of the suspension (flexible rotary cable principle) from 
its normal. The upper and lower ends of the fiber suspensions are turned 
through angles y and by the normal torque component. 

There are five generalized coordinates - \p, cf>, Yr fi/ and 0 
represented by q^. and dissipative generalized forces due to damping. 
According to Lagrange's equation, namely. 


d 

dt 



3a 

8qr 


Generalized force 


(B116) 


five equations of motion are obtained which can be solved, as was done in this 
paper, for the output with perturbation effects. 
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DERIVATION OF DAMPING COEFFICIENT FROM A DECAY CURVE USING 
THE AXIS OF SYMMETRY AS THE REFERENCE DATUM 

In this appendix, the equations for damping coefficient c and rigidity K 
are derived. To apply these equations, only two peak torsional amplitudes 
(measured from the axis of symmetry of the decay curve) and the number of half- 
periods between them must be known. 

The equation for an unperturbed, freely oscillating torsion pendulum is 


J0 + £0 + K0 = 0 


(Cl) 


where J is the polar torsional moment of inertia, £ is the damping constant, 
K is torsional rigidity, and 0, 0, and 0 are the angular displacement, 

velocity, and acceleration of the pendulum. Let 


and 


b 


Jt 

2J 


0 ) 


K 

J 


(C2) 


(C3) 


Then equation (Cl) becomes 


0 + 2b0 + co^ 


= 0 


(C4) 


and the solution of equation (C4) is 


exp^-b + \j^ 


2 2 

W yt + c^ exp 


(-b - \|b^ - (x/)t 


The initial conditions at t = 0 are 


= 0 


(C5) 


(C6) 
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With these initial conditions, the constants C]^ and C 2 can be obtained, so 
that equation (C5) becomes 


0o 

0 = — exp (-bt) 


+ 1 - 


1 + 


L2 2 

b - 0) 


exp ( \| b^ - 0)^ t) 


..2 


j exp(- 


- \lb^ - t 


/b'- - OJ' 

of interest is the underdamped case when 


(C7) 


b < (jQ 


= i \f(j 


io)' = i - b^ 


(C8) 


where oj' is the damped angular velocity. Equation (C7) becomes 


9 = 0Q exp(-bt)^cos OJ't + ~ sin 00 ' t^ 


(C9) 


When t equals the period of oscillation T, 


co'T = 2 tt 


(CIO) 


and from equations (C2) , (C3) , (C8) , and (CIO) 


Let 


rn — . 


2TT 


i — 



1/2 


^ - ( 

' £ 



J \ 

.2jj_ 


b _ 


£ 


03 

2(KJ) 

equation 

(C12) 

in equation 


(Cll) 


(Critical damping occurs when c = 1.0) (C12) 


LO 


p 1/2 

(1 - c^) w 


(C13) 


and from equations (Cll) and (C2) 


K = 


4TT^J 


+ c2q)2j 


(C14) 
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The phase angle 6 is 


^ , -1 
0 = tan — 
b 


(C15) 


Equations (C8) , (C9) , (C12) , (C13) , and (CIS) give 


0Q exp(-c03t) 
® ~ - 1 /o 


JT7T 

(1 - C^) 


t + 6) 


(C16) 


The phase angle does not affect the period. 

A typically damped output, on an expanded time base, would look like fig- 
ure 11. To obtain the value of c from this output, initial conditions must be 
taken at a time t when 0=0. Hence, at any arbitrary peak displacement, a 
time reference t = 0 is taken to obtain from equation (C16) , 


sin 6 


(C17) 


(1 - c^) 


At t = Vn/uj' , where n is the number of half-periods between 0^^ and 0^^, 


(1 - c^) 


/-CLOTTn) , ~ 

expl I sin 0 


(C18) 


From equations (C17) and (C18) , 


0 . / 

1 / COJ 

^ = exp^— tn 
n ' 


(n = 0,1,2, . . 


(C19) 


In (0i/0j,) 


(C20) 


Substituting equation (C13) results in 


1 + a-^ 


(C21) 
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For obtaining c, the damping coefficient of the material, K, its 
rigidity, and G, its shear modulus, only two peak points and the number of 
half-periods separating them are needed. Since the expression used for obtain- 
ing c (eq. (C19) ) is the decay envelope of the peaks, if two alternate peaks 
are used, the algebraic sign of one of them must be reversed. Either way, the 
axis 0=0 must be accurately known; this is a great disadvantage. In 
appendix D, a three-successive-peak method is derived which does not have this 
drawback . 
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THREE-SUCCESSIVE-PEAK METHOD OF OBTAINING DAMPING COEFFICIENT 
FROM A DECAY CURVE USING ANY DATUM AS A REFERENCE 

In appendic C, equations were derived for obtaining damping coefficient c 
with two peak amplitudes and the number of half-periods between them known. 
However, with this method the peak values must be measured from 0=0, the 
axis of symmetry of the decay curve. There are numerous other ways of obtain- 
ing c, some employing alternate peak measurements and others employing succes- 
sive peak measurements. It is felt that three successive peak measurements are 
best from the point of view of accuracy. One of the reasons for this conclusion 
can be seen from the calibrate curve in figure 13. If three successive points 
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Figure 13.- Calibrate curve showing galvanometer output as a 
function of pendulum window angle input. 
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are taken, the effects of nonsymmetries, variations in sensitivity, and non- 
linearity are either minimized or eliminated. In this appendix, the three- 
successive-peak method is derived; its only requirement is that the reference 
axis from which the three successive peak values are measured be parallel to 
0=0. This is especially advantageous from the point of view of computer 
readout of data. 

From equation (CIS) , 0 for successive peaks is 


(1 - c^) 


1/2 


caj2TTn\ . r n ( CL02im\ 

exp - ; — I sin 0=0 exp - : 

' W / o ^ V w' 




(n = 0,1,2, . . .) (Dl) 


1/2 

where sin 6 = 0)'/w and u)' = (1 - c^) ca. Also, 


n+1 


= 0 


exp 


CO) 

co' 


21T (n + 1) 


) 


(D2) 


Consider the ratios of the general terms in equations (Dl) and (D2) : 

^ COJ \ 

1 - — r 2TTn 


exp 


n+1 




COJ 


-o 2u(n + 1 


) 


= exp(^ 2^) = 


Constant 


(D3) 


Now consider a translation of the 0,t coordinate system to the 9',t coordi- 
nate system, as shown in the following sketch: 



d 
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For this coordinate system 


e; = 9n - d 


(n = 0,1,2,...) (D4) 


Substituting equation (D4) in equations (Dl) and (D2) results in 


0 ' + d = 0 expf- 2Trn^ 
n o 03' / 


(D5) 


®n+l ^ 


d = 0^ exp 27T (n + l)j 


(D6) 


The ratio of equations (D5) and (D6) in the 0',t system is the same as in the 
0,t system, that is. 


+ d 


exp 


( CO) ^ \ 


e + d 

n+1 


C/i) 

exp - — - 2ir(n + 1) 

L 


/co) „ , 

= expl 2TT j = Constant 


(D7) 


This equation can be solved for c to result in the two-point method 
derived in appendix C. However it is dependent on d. To avoid measuring d 
in equation (D7) , consider the next successive peak 


^n+2 + = 60 


exp [j 2 tt (n + 2)j 


^n+1 ^ ^ /co) ^ 

— = exp — r 2 tt 

>n*2 * « 


) 


(D8) 


From equations (D7) and (D8) 


/ COJ 


+ d = exp(^ 2 tt + d) 


(D9) 


n+1 


+ d = exp(^ 27r)(e;_^2 + 


(DIO) 
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Subtract 9^+]^ + from both sides of equation (D9) and 9^+2 from both 

sides of equation (DIO) to obtain 



(Dll) 


(D12) 


Taking the ratio of equations (Dll) and (D12) and substituting equation (D8) 
result in 


9n - 9ntl 


®;ti 




Qn+1 + ^ 

0' ^ + d 
n+2 



For n = 0 , 


e; - 0' 


exp 




Taking the logarithm of both sides of equation (D14) results in 


cgj 

CO' 



a 


(D13) 


(D14) 


(D15) 


This is similar to equation (C20) . Note that this method requires three 
points 0Q, 0j, and 0^ measured from any reference axis parallel to 0=0. 

Since this method is independent of the position of 9' = 0, it is very desir- 
able. No restriction was placed on the 9 values; thus, they are the measured 
values including their sign as measured from the reference datum. For these 
reasons, this method is used for obtaining a and, from equation (C21), c, the 
dynamic damping. Note that the above development for three successive peaks may 
be readily extended to n successive peaks. For alternate peaks, however, the 
sign of the alternate peaks has to be reversed so that they lie on the exponen- 
tial curve. Hence, for three alternate peaks. 


CO) 

U)' 


In 


' e; - Cl \ 

0Atl + 9At2/ 


IT 


(D16) 


Again, the algebraic values of 0' including their sign are measured from the 
reference datum. 
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